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ABSTRACT 

We present near-infrared K-band spectroscopy of 21 elliptical or cD Brightest Cluster 
Galaxies (BCGs), for which we have measured the strength of the 2.293 /im CO stellar 
absorption feature. We find that the strength of this feature is remarkably uniform 
among these galaxies, with a smaller scatter in equivalent width than for the normal 
elliptical population in the field or clusters. The scatter for BCGs is 0.156 nm, com- 
pared with 0.240 nm for Coma cluster ellipticals, 0.337 nm for ellipticals from a variety 
of other clusters, and 0.422 nm for field ellipticals. We interpret this homogeneity as 
being due to a greater age, or more uniform history, of star formation in BCGs than 
in other ellipticals; only a small fraction of the scatter can be due to metallicity vari- 
ations, even in the BCGs. Notwithstanding the small scatter, correlations are found 
between CO strength and various galaxy properties, including R-band absolute mag- 
nitude, which could improve the precision of these galaxies as distance indicators in 
measurements of cosmological parameters and velocity flows. 
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1 INTRODUCTION 

Elliptical or cD Brightest Cluster Galaxies (BCGs) have long 
been used as cosmological probes, in attempts to determine 
parameters such as Ho and qo (Sandage 1972; Sandage & 
Hardy 1973; Hoessel, Gunn & Thuan 1980; Sandage & Tam- 
mann 1990). These studies make use of the ease of selection 
of such objects at cosmological distances, which is a result 
of their bright optical luminosities and privileged locations 
defined by other galaxies and by the centroids of cluster 
X-ray emission. They also appeal to another, more surpris- 
ing property of BCGs: they appe ar to b e excellent standard 
candles, with Sandage & Hardy ( 1973 ) finding a scatter of 
only 0.28 mag in visual ab solute magnitude, after correction 
for cluster Bautz-Morgan (197C) and richness class. Lauer & 
Postman (1994), in a study using more modern CCD tech- 
niques, found a scatter of 0.33 mag in R-band absolute mag- 
nitudes of BCGs. They were able to reduce this scatter to 
0.25 mag by using the relation of absol ute m agnitude with 
the 'Structure Parameter' 01 of Hoessel (1980), which is de- 
fined as 



d log L„ 



Q as being a dimensionless parametrization of galaxy size, 
and empirically it is found to increase with BCG luminos- 
ity. Indeed he found a scatter in a-corrected BCG absolute 
magnitudes of jus t 0.21 mag, even lower than that of Lauer 
& Postman ( |l99^ ). 

This uniformity in BCG luminosities remains something 
of a mystery. There is indeed good reason for expecting a 
wider disparity in BC G pr operties than in t hose of other 
galaxy types. Hoessel (1980) and Lauer (L988) showed that 



dlog r 



where Vm is a metric radius of 9.6 kpc and Lm is the lu- 



minosity within that metric radius. Hoessel ( 1980 ) explains 



~30% of BCGs have multiple nuclei, a frequency which is 
greater than would be predicted by chance superpositions, 
and which indicates that galactic cannibalism is very com- 
mon in BCGs at the present epoch. The continued accre- 
tion of cluster members would seem to ensure that BCGs 
as a population should have diverse properties, and it would 
seem probable that this should be reflected both in their 
total luminosity and in their composition. One motivation 
for the present paper is to test whether BCGs are more ho- 
mogeneous or more diverse than the general population of 
elliptical galaxies. 

In this paper, we look at the stellar populations of 
BCGs, by measuring the strength of the 2.293^m CO ab- 
sorption feature for 21 BCGs, and comparing the distribu- 
tion of values with similar measurements for a large sample 
of field, group and cluster ellipticals (Mobasher & James 
1996; James & Mobasher 1999; Mobasher & James 2000). 
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CO strength contains information on recency of star forma- 
tion, since it is very strong in supergiants (present lO'^-lO* 
years after a burst of star formation) , strong in the cool AGB 
stars which contribute significantly to the near-IR light after 
10**-10^ years (Renzini & Buzzoni 1986; Oliva et al. 1995), 
and somewhat weaker in older populations. It also displays 
some metallicity depen dence, being weak in very low metal- 
licity globular clusters (Origlia et al. 1997). This dependence 



was quantified by Doyon, Joseph & Wrigh t (1994), and fur- 
ther studied in Mobasher & James ( |l999|) . 



Such studies of BCGs are particularly significant consid- 
ering the apparent large-scale velocity flow found by Lauer 
& Postman (1994). Using a sample of 119 BCGs out to a 
redshift of 15000 kms~^, they found the restframe defined 
by the galaxies to differ from that of the Cosmic Microwave 
Background by almost 700 kms^^. This result has been in- 
terpreted as evidence for a cosinological streaming fiow, but 
an alternative explanation would be that BCG properties 
vary systematically around the sky, for example due to stel- 
lar population changes from galaxy to galaxy. This provides 
a further motivation for the present study. 

The organisation of this paper is as follows. Section 2 
describes the selection of target galaxies, the observations, 
and the data reduction. Section 3 contains the main re- 
sults, including a comparison of CO absorption strengths 
of BCGs and other elliptical galaxies, and correlations of 
CO strengths with other galaxy parameters. Section 4 sum- 
marises the main conclusions. 



2 SAMPLE SELECTION, OBSERVATIONS 
AND DATA REDUCTION 

The galaxy s ample was selected from the BCG list of Lauer 
& Postman ( 1994 ). All have measured recession velocities 
less than 15,000 kms~^, with R band photometry presented 



in Lauer & Postman ( 1994 ) . The observations presented here 
were carried out using the United Kingdom Infrared Tele- 
scope (UKIRT) during the 4 nights of 21-24 February 1999. 
The instrument used was the long-slit near-IR spectrometer 
CGS4, with the 40 line mm~^ grating and the long-focal- 
length (300 mm) camera. The 4-pixel-wide slit was chosen, 
corresponding to a projected width on the sky of 2.4 arcsec. 
Working in 1st order at a central wavelength of 2.2 ^m, this 
gave coverage of the entire K window. The CO absorption 
feature, required for this study, extends from 2.293 /im (rest 
frame) into the K-band atmospheric cut-off. The principal 
uncertainty in determining the absorption depth comes from 
estimating the level and slope of the continuum shortward 
of this absorption which requires wavelength coverage down 
to at least 2.2 fj,m and preferably to shorter wavelengths. 
There are many regions of the continuum free from lines 
even at this relatively low resolution. The effective resolu- 
tion, including the degradation caused by the wide slit, is 
about 230. 

For each observation, the galaxy was centred on the slit 
by maximising the IR signal, using an automatic peak-up 
facility. Total on-chip integration times of 12 minutes were 
used for the brightest and most centrally-concentrated el- 
lipticals while an integration time of 24 minutes was more 
typically required. During this time, the galaxy was slid up 
and down the slit at one minute intervals by 22 arcsec, giv- 



ing two offset spectra which were subtracted to remove most 
of the sky emission. Moreover, the array was moved by 1 
pixel in the spectral direction between integrations to en- 
able bad pixel replacement in the final spectra. Stars of 
spectral types A0-A6, suitable for monitoring telluric ab- 
sorption, were observed in the same way before and after 
each galaxy, with airmasses matching those of the galaxy 
observations as closely as possible. Flat fields and argon arc 
spectra were taken using the CGS4 calibration lamps. A to- 
tal of 21 brightest cluster galaxies was observed. 

The data reduction was performed using the FIGARO 
package in the STARLINK environment. The spectra were 
flatfielded and polynomials fitted to estimate and remove the 
sky background. These spectra were then shifted to the rest 
frame of the galaxy, using redshifts from Lauer & Postman 
(1994). The atmospheric transmissions were corrected by di- 
viding each spectrum with the spectrum of the star observed 
closely in time to the galaxy, and at a similar airmass. The 
resulting spectrum was converted into a normalised, recti- 
fied spectrum by fitting a power-law to featureless sections 
of the continuum and dividing the whole spectrum by this 
power-law, extrapolated over the full wavelength range. Two 
rectified spectra are shown in Fig. 1. The apparent emission 
features at 2.14 fim and 2.10 fim are artefacts caused by ab- 
sorptions in the A stars used for atmospheric transmission 
correction, and appear in different positions because of the 
restframe corrections. 

To measure the depth of the CO absorption feature, 
the procedure outlined in James & Mobasher (1999) is used. 
The restframe, rectified spectra were rebinned to a common 
wavelength range and number of pixels, to avoid round- 
ing errors in the effective wavelength range sampled by a 
given number of pixels. The CO strength for each spec- 
trum was deter mined using the method of Puxley, Doyon 
& Ward (1997). They advocate the use of an equivalent 
width, CO Ew, which is determined within the CO absorp- 
tion feature between rest-frame wavelengths of 2.293 /im 
and 2. 320 /x m. This wavelength range was found by Puxley 
et al. (1997) to give maximum sensitivity to stellar popula- 
tion variations, and can be used for galaxies with recession 
velocities of up to ~ 18000 kms~^ before the spectral region 
of interest shifts out of the usable K window. This is not the 
case for the CO index COs 



jsp used by Doyon et al. (1994) 



which extends over a restframe wavelength range of 2.320- 
2.400 and would have been affected by large and uncer- 
tain telluric absorption and emission for the highest redshift 
galaxies in the present sample. Thus we only present CO ew 
values in this paper. 

A furth er advantage of the CO ew definition of Puxley 
et al. (1997) is that it is almost completely unaffected by 
velocity dispersion effects, due to the wide range of wave- 
leng th ov er which the absorption is measured. Puxley et 
al. (1997) find the velocity dispersion corrections to be in- 
significant, which we confirmed by smoothing low- velocity- 
dispersion galaxy spectra to an effective velocity dispersion 
of 500 kms~^. The resulting change in COew was ~0.25%, 
very much smaller than the random errors. 

The errors on the COew values include three compo- 
nents. The first was calculated from the standard deviation 
in the fitted continuum points, on the assumption that the 
noise level remains constant through the CO absorption, giv- 
ing an error on both the continuum level and on the mean 
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Figure 1. The rectified spectra of the BCGs in Abell 1060 (top) 
and Abell 1016 (bottom), showing the wavelength ranges used to 
define the continuum level and the CO Equivalent Width. These 
spectra illustrate the lowest and highest signal-to-noise spectra 
used for the present paper. 



level in the CO absorption, which were added in quadrature. 
The second error component comes from the formal error 
provided by the continuum fitting procedure. This proce- 
dure could leave a residual tilt or curvature in the spectrum, 
and the formal error was used to quantify this contribution. 
The final component was an estimate of the error induced by 
redshift and wavelength calibration uncertainties. All three 
errors were of similar sizes, with only the first varying from 
spectrum to spectrum, as a result of signal-to-noise varia- 
tions (see Fig. 1), and all three were added in quadrature to 
give the value quoted in Table 1. 





Figure 2. Histograms of COew distributions BCGs (hashed re- 
gions), overlaid on field and group galaxies (a), cluster galaxies 
(b), and Coma cluster galaxies (c). 



3 RESULTS 

The equivalent widths of CO absorption features for the 
sample of BCGs observed in this study are presen ted in 
Table 1. The data included in this table are Abell (1950 



catalogue numbers (column 1), BCG names (column 2), 
COew values with l-cr errors (column 3), recession veloc- 
ity in kms~^ (column 4), absolute R-band magnitude corre- 
sponding to the metric luminosity Lm (column 5), structure 
parameter (a) (column 6) and the magnitude residual rela- 
tive to the best-fit L,n-a relation (column 7), (columns 4-7 



are all taken from Lauer & Postman (1994), who assumed a 
Hubble constant of 80 kms^^Mpc"^). Columns 8 and 9 con- 
tain velocity dispersions and Mg2 metallicity indices, where 
available, from Faber et al. (1989). 

We find the mean COew value for the 21 BCGs 
(3.35±0.03) to be effectively id entical to that of the Corn a 
cluster ellipticals (3.37±0.04) ( |Mobasher fc James 1999|) , 



and to that of the 31 ellipticals f rom a range of clusters 
discussed by James & Mobasher ( |l999| ) (3.29±0.06). The 
cluster and BCG distributions lie between the distributions 
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Table 1. Photometric and spectroscopic parameters for 21 brightest cluster galaxies 
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of 'isolated' a nd 'gr oup' field ellipticals discussed by James 
& Mobasher (1999) and shown in Fig. 2a. The major dif- 
ference between the BCG COew values and those of other 
ellipticals is the remarkably small range in the former: the 
standard deviation for BCGs is 0.156, compared to 0.240 for 
Coma ellipticals, 0.337 for general cluster galaxies, and 0.422 
for cluster plus field ellipticals. Indeed, the scatter in BCG 
CO absorption strengths is that predicted from the error es- 
timates on the individual COew values, and so the intrinsic 
scatter may be much smaller still. Given the small number 
of BCG galaxies, a Kolmogorov-Smirnov test cannot distin- 
guish between the distributions of BCG and cluster or Coma 
galaxies in Figs. 2b and 2c, but there is less than 10% chance 
that the BCGs are drawn from the same parent population 
as all the non-BCG ellipticals, and less than 1% chance that 
they are from the same population as field ellipticals (Fig. 
2a). 

The BCGs are drawn from a much narrower region of 
the galaxy luminosity function than are the comparison sam- 
ples in Fig. 2, which could affect the interpretation of this 
resuh. The 21 BCGs have a range in Mr of -22.0 to -23.1, 
little more than a magnitude. R-band photometry is not 
available for all the comparison galaxies, but good estimates 
can be made from published optical and near-IR photome- 
try, leading to an estimated range of M_r of -19.4 to -22.8 
for the Coma cluster ellipticals, and -19.8 to -22.5 for the 
field a nd cluster ellipticals discussed by James and Mobasher 
(1999). We investigated whether the differences in COew 



scatter shown in Fig. 2 result from these differences in lu- 
minosity range by regressing CO ew on absolute magnitude, 
and studying the distributions of COew residuals about the 
best-fit lines. The distributions of these residuals are shown 
in Fig. 3. The dashed, diagonally shaded columns represent 
the residuals for the BCGs; the thick, dotted columns are 
those for the Coma cluster galaxies; and the solid lines rep- 
resent the residuals for the field, group and cluster galaxies 




RESIDUAL ACO„ 



Figure 3. Histograms of residuals of CO^w when regressed on 
absolute magnitude, for BCGs (hashed regions with dashed hnes), 
overlaid on residuals for Coma ellipticals (thick dotted lines) and 
general field and cluster ellipticals (solid lines). 



from James and Mobasher (1999). The standard deviations 
of the COew residuals are 0.133 nm for the BCGs, 0.221 nm 
for the Coma ellipticals, and 0.419 nm for the cluster plus 
field ellipticals. This reinforces the conclusion from Fig. 2 
that the BCGs have substantially more homogeneous CO 
strengths than the other elliptical galaxies studied, and this 
result does not appear to be a selection effect caused by the 
small luminosity range of the BCGs. 

This uniformity in COew values is the main result of 
this paper, and it is important to consider what it implies 
in terms of differences between BCGs and other ellipticals. 
Both high metallicity and recency of star formation are ex- 
pected to increase COew values. The effect of metallicity on 
COew values can be estimated for the galaxies with mea- 
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sured Mg2 indi ces, u sing the following method. From Fig. 
37 of Worthey (|l994|), a change in [Fe/H] from -0.25 to 0.00 



changes the Mg2 index from 0.216 to 0.258, and the change 
is approximately linear over the modelled range. Thus we 
infer a relation of the form 

SMg2 = 0.168 S[Fe/H]. 



Doyon et al. (1994) find a relation between [Fe/H] and their 
CO index CO^p, 

SCOsp = 0.11 5[Fe/H], 



and from the definitions in Puxley et al. ( 1997 ) it is straight 



forward to convert from the index CO sp to CO ew- Then, 
the measured scatter in Mg2 index of 0.029 for the BCGs 
should cause a scatter of 0.060 nm in COew, 38% of the ob- 
served scatter. For Coma ellipticals, the measured Mg2 scat- 
ter is 0.024, equivalent to a scatter of 0.049 nm in COew, 
20% of that observed, and for the field and cluster sample, 
the Mg2 scatter is 0.030, and the predicted CObvi' scatter 
0.062 nm, 15% of that observed. Note also that the scat- 
ters in Mg2 values are very similar in the three subsam- 
ples, whereas they have very different COew distributions. 
Thus, we conclude that metallicity differences have little ef- 
fect on the measured CO ew values for the elliptical galaxies 
studied here, and propose that star formation history is the 
dominant factor causing the larger scatter for non-BCG el- 
lipticals. If so, the differences in the distributions of COew, 
shown in Fig. 2 would be the result of wider variations in 
star formation history for general field and cluster ellipticals 
than for the BCGs. This indicates that BCGs formed their 
stars very early; if there has been more recent star forma- 
tion in these galaxies then the rate of star formation as a 
function of epoch must have been very uniform from galaxy 
to galaxy. 

Given the narrow range in BCG COew values, it is un- 
realistic to expect very strong correlations with other BCG 
parameters. Nevertheless, Fig. 4 does show a good correla- 
tion with absolute R-band magnitude in a 10 kpc metric 
aperture, M_b, with a correlation coefficient of 0.51 and a 
probability of 98.4% that this represents a true correlation 
(i.e. 1.6% probability that it could arise by chance). This 
is significant enough to be useful as a distance indicator: 
the scatter in M_b for the 21 galaxies observed is 0.326 mag, 
which reduces to 0.280 mag when the Mr values are cor- 
rected for the COew effect. The slope of the regression line 
of M_R on COew is somewhat smaller than that for the 
trend in Mk (total K-band absolute magnitude) vs C Oew 
for Coma cluster ellipticals (Mobasher & James 1999), at - 
1.1±0.5 mag/nm c.f. -1.6±0.7 mag/nm for the Coma galax- 
ies. However, this difference is not statistically significant 
(~0.6a). It is not possible to determine whether the abso- 
lute magnitude-CO ew relations are consistent for the vari- 
ous samples because of the lack of homogeneous photometry, 
and the consequent need for large and uncertain colour and 
aperture corrections. 

Similarly, there is a strong correlation between COew 
and residuals (dM ^) about the relation of M_r with struc- 
ture parameter a (Hoessel 1980) (Fig. 5), in the sense that 
galaxies with high COew tend to be bright relative to 
the mean relation (correlation coefficient 0.60, significance 
99.6%). The residuals (dM^) are reduced from 0.243 mag to 
0.195 mag by correcting for the trend with COew shown 




3.2 3.4 
CO Equivalent Width 



3.6 



Figure 4. Absolute R-band magnitude versus CO Equivalent 
Width for the 21 Brightest Cluster Galaxies. 




3.2 3.4 
CO Equivalent Width 

Figure 5. Structure parameter residuals versus CO Equivalent 
Width for the 21 Brightest Cluster Galaxies. 



in Fig. 5. There is no correlation between COew and the 
structure parameter a itself. 

Given the trends found in Figs. 4 & 5, it is instructive 
to explore if these effects could cause the putative streaming 
flow signal detected by Lauer & Postman (1994) using the 
full sample of BCGs. However, we find no significant cor- 
relation between COew and direction o n the sky (Fig. 6). 
This implies that the Lauer & Postman (1994) apparent de- 



tection of a bulk flow was not an artefact of differing stellar 
populations between sample galaxies, although we have of 
course only looked at a small fraction (18%) of their sample. 

The trends in flgures 4 & 5 are both in the sense that 
brighter galaxies have higher indices: that in Fig. 4 could be 
a consequence of a metallicity-absolute magnitude relation, 
and there is indeed evidence of a weak correlation of CO ew 
with metallicity (Fig. 7). The correlation coefficient here is 
0.52, but the relation has only 80% significance due to only 
8 BCGs having tabulated Mg2 values. Fig. 7 also shows the 
t rend in COew with meta. llicitv for 31 Coma cluster galaxies 
(Mobasher & James 199!:). The BCGs clearly lie at higher 
mean metallicity than do the Coma cluster ellipticals (mean 
Mg2 values 0.328 for the BCGs and 0.294 for the Coma el- 
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Cos(9) 



3.5 

CO Equivalent Width 



Figure 6. CO^w versus the cosine of the angle betw een galaxy 
direction and the apex of the Lauer & Postman(1994) streaming 
motion, for the 21 Brightest Cluster Galaxies. 



Figure 7. Mg2 index versus CO Equivalent Width for the 
Brightest Cluster Galaxies (circles) and Coma cluster ellipticals 
(crosses). 



lipticals). Using the relation between Mg2 and [Fe/H] from 
Worthey ( 1994 ), and that between [F e/H] and CO absorp- 
tion strength found by Doyon et al. (1994), this predicts a 
difference in mean CO ew of 0.05 nm, compared to the ob- 
served difference of 0.04±0.07 nm in the same sense. The 
small variation found here confirms the weakness of the de- 
pendence of CO EW on metallicity, as found by Doyon et 
al. (1994), at least at the high metallicity values typical of 
centres of bright galaxies, and also confirms our earlier con- 
clusion that star formation history is the dominant effect in 
determining CO ew strength. This is further confirmed by 
the wide s pread in the CO ew valu es of isolated and group 
ellipticals ( James fc Mobasher 199£| ), with no corresponding 
change in their metallicity. 

Finally, we find a weak correlation of CO ew with veloc- 
ity dispersion for 14 galaxies with data in Table 1 (Fig. 8). 
This corresponds to a correlation coefficient of 0.294, and 
is significant at the 69% level. The slope and correlation 
coefficient are the same as was found for 31 Coma cluster 
ellipticals, also plotted in Fig. 8, but the mean correlation 
for the BCGs is again offset, to higher velocity dispersion at 
a given CO ew- The mean CO ew is almost identical for the 
two samples, whilst the BCGs have a much higher average 
velocity dispersion, and hence mass, as expected. 



4 CONCLUSIONS 

We find that BCGs are much more homogeneous in evolved 
red stellar content than ellipticals overall, and BCGs are 
somewhat more homogeneous than Coma cluster ellipticals. 
The measured scatter in the CO ew indices for BCGs is 
comparable to the measurement errors. We interpret this as 
implying a more uniform and probably earlier star formation 
history for BCGs than for normal ellipticals. Metallicity does 
not appear to be the controlling parameter of CO absorption 
strength for elliptical galaxies. 

Absolute magnitudes, and magnitude residu als re lative 
to the structure parameter relation of Hoessel (198C) cor- 
relate well with CO absorption depth. This may imply the 
presence of an additional intermediate-age population, or 




3.5 

CO Equivalent Width 



Figure 8. Log of velocity dispersion (kms~^) versus CO Equiva- 
lent Width for the Brightest Cluster Galaxies (circles) and Coma 
cluster ellipticals (crosses). 



a higher metallicity population, in the galaxies which are 
over-luminous relative to the mean relation. However, this 
effect is very small and such populations, if present, must be 
weaker than in most field and cluster ellipticals, given the 
high degree of homogeneity in BCG CO ew values. 

It may be possible to use these correlations to define 
higher precision distance indicators for BCGs, which could 
be used for galaxies with redshifts up to 18,000 kms~^, as is 
indicated by the reduced scatter in R-band absolute magni- 
tude, and the reduced scatter about the structure parameter 
relation after correction for the correlation with CO ew, dis- 
cussed in section 3. 

Recent ROSAT observations reveal that a significant 
number of the Lauer & Postman (1994) galaxies do not lie 



at the X-ray centroids of their clusters (Paul Lynam, private 
communication), and there may be better candidates for 
the dominant central cluster galaxy. Thus our conclusions 
may refer more generally to bright galaxies towards cluster 
centres than to individual BCGs inhabiting the very centre 
of the cluster potential. 
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